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Point defects are ubiquitous in two-dimensional crystals and play a fundamental role in deter- 
mining their mechanical and thermodynamical properties. When crystals are formed on a curved 
background, finite-length grain boundaries (scars) are generally needed to stabilize the crystal. We 
provide a continuum elasticity analysis of defect dynamics in curved crystals. By exploiting the fact 
that any point defect can be obtained as an appropriate combination of disclinations, we provide 
an analytical determination of the elastic spring constants of dislocations within scars and compare 
them with existing experimental measurements from optical microscopy. We further show that va- 
cancies and interstitials, which are stable defects in flat crystals, are generally unstable in curved 
geometries. This observation explains why vacancies or interstitials are never found in equilibrium 
spherical crystals. We finish with some further implications for experiments and future theoretical 
work. 

PACS numbers: 82.70.Dd,61.72.Bb,61.72.Ji,61.72.Mm 



I. INTRODUCTION 

The rich physics of the ordering of matter on planar 
surfaces takes on a new complexion when the ordering 
occurs on a curved two-dimensional manifold. Gaussian 
curvature, for example, favors the appearance of topo- 
logical defects that are energetically prohibitive in planar 
systems. This has been demonstrated in the case of suf- 
ficiently large spherical crystals [H, H, S 0> H> > toroidal 
hexatics Q, and both crystals and hexatics draped over 
a Gaussian bump [8|, l9( . 

For the simplest case of crystalline order on the con- 
stant curvature two-sphere (the surface of a solid ball 
in R 3 ) the key new feature is the appearance of scars 
(Fig. [I]), linear strings of dislocations around a cen- 
tral disclination that freely terminate inside the crystal, 
for crystals with radius above a microscopic-potential- 
dependent critical radius [l], 0, H[ ■ Scars have been ob- 
served experimentally in systems of colloidal beads self- 
adsorbed on spherical water droplets in an oil emulsion 
[To| . The imaging technique (conventional microscopy) 
in these experiments only allowed spherical caps covering 
10% — 20% of the full sphere to be imaged. Recently the 
use of fluorescently labeled colloidal particles and laser 
scanning confocal microscopy allowed the imaging of 50% 
of the sphere. In this way the global spatial distribution 
of scars was also measured [111 ]. 

Recent experiments fl2l [l3| have investigated the dy- 
namics of defects by directly visualizing colloidal parti- 
cles absorbed on spherical oil- water interfaces. It was 
shown that dislocation glide within the scars (see Fig. [2]) 
could be described very accurately by a harmonic poten- 
tial binding the dislocation to the scar and an empirical 
Peierls potential that models the underlying crystalline 
lattice. The spring constants of the harmonic potentials, 
the elastic stiffness of the dislocation, were obtained from 
fits to the experimental results. In this paper, we show 



that continuum elasticity theory [l], 0] can be used to 
provide explicit first-principles predictions for the elastic 
stiffness. 

Defects such as vacancies and interstitials are quite 
common in two-dimensional crystals [l4]. It has been 
shown that, quite generally, the presence of vacancies and 
interstitials significantly reduces the crystal's strength as 




FIG. 1: (Color online) A light microscope image, reproduced 
from [12|. of a water droplet with an 85 /jm diameter and 1.9 
fj,m mean particle spacing (R/a « 22). Fivefold (+1) disclina- 
tions are colored red, Sevenfold (—1) disclinations are colored 
yellow, and tightly bound five to seven pairs represent dislo- 
cations. The three dislocations whose dynamics are analyzed 
in this paper are displayed within the blue box. 
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a result of stress enhancement effects [HI . Jain and Nel- 
son (l6j performed an extensive investigation of inter- 
stitials and vacancies in two-dimensional planar crystals 
and identified three different intcrstitials and vacancies, 
depending on their symmetry, as the prevalent struc- 
tures. Subsequent experiments [171 ] confirmed the stabil- 
ity of these defects and studied their dynamics. Very re- 
cently, Brownian dynamics simulations [181 ] have revealed 
a complex kinetics with a variety of modes that allow 
defects to glide and rotate. Rather interestingly, vacan- 
cies and interstitials have not been observed either ex- 
perimentally or in numerical simulations [1 91 ] in spherical 
crystals. In this paper, we provide a study of the stability 
of vacancies and interstitials in curved two-dimensional 
crystals. Our analysis uses continuum theory, and there- 
fore the results are directly applicable to other systems 
such as, for example, the analysis of vacancies and their 
relation to failure stress, which has recently investigated 
in straight carbon nanotubes [20| , which provides another 
example for curved crystals. 

The paper is organized as follows: In Sec. |TT] we de- 
scribe the dynamics of scar defects and the continuum 
elasticity theory of defect interactions. The theoretical 
results so obtained are compared to the experiment in 
Sec. IIII1 In Sec. [TV] we study the instability of intersti- 
tials and vacancies in curved crystals with the continuum 
elastic model. The Thomson problem java applet used 
for this analysis is described in the Appendix. 
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FIG. 3: (Color online) Spring model of dislocation binding. 
The red dot (D) represents the central disclination and the 
blue sticks (i = 1,2,...) represent successive dislocations em- 
anating from the central disclination. The spring constants 
fci, &2,... represent the binding of each dislocation to its parent 
scar. 



II. DYNAMICS OF SCAR DEFECTS 

A. Empirical description of scar dynamics 

The scar dynamics is obtained from light microscopy 
as discussed in Ref. 12]. A typical snapshot of a config- 
uration is shown in Fig. [1] and its evolution as a function 
of time is shown in Fig. [2J showing some dislocations 
within the scar gliding at different times. In Ref. [Hj], it 
was shown that these data are well described by a model 
where each dislocation within the scar is pinned by a har- 
monic potential with spring constant ki (here i labels the 
position of the dislocation within the scar) as shown in 
Fig. [3] — that is, 
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FIG. 2: (Color online) The dynamics of the three dislocations 
within the blue box of Fig. [1] taken from Ref. [l^]- Motion 
consists of local Brownian fluctuation about the equilibrium 
position together with larger-scale glide. 



where Sj is the geodesic displacement of the zth dislo- 
cation on the surface of a sphere and the last term is 
the Peierls potential [2l[, which models the underlying 
crystalline structure of the lattice. Values for the exper- 
imentally determined spring constants were determined 
m Ref. :ia. We now provide the details leading to an 
explicit evaluation for the elastic stiffness. 



B. Continuum elasticity of scars in curved 
backgrounds 

We first present a discussion of the elasticity of scars. 
Point topological defects can be parametrized by discli- 
nations. We therefore introduce a disclination density 
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^Qi(5(x - Xi) , 
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where qt is the disclination charge (qi = +1 for 5's and —1 
for 7's). The elastic energy of an arbitrary disclination 
density has been discussed extensively 0, [H, Q and is 
given by 

E = j f j <f^d 2 y^)^(y)[K^)-Q^)] -L 

- ■' -> X y 
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x[K(y)-Q(y)]+NE c 



(3) 



where K(x) is the Gaussian curvature of the background 
with metric g(x), Y is the two-dimensional Young's mod- 
ulus, and E c is the disclination core energy. Both Y and 
E c depend on the microscopic particle potential. 

The free energy of Eq. ([3]) for a spherical crystal with 
the disclination density of Eq. © is then [T], @, y] 
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where 
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and (3 is the angular geodesic length between points 
(9i,(f>i) and (9j,<pj): 



cos j3 = cos 8i cos 9j + sin 6i sin 9j cos( 



■ (6) 



The previous energy is a function of disclinations only. 
It is convenient to introduce dislocations explicitly. A 
dislocation can be regarded as a tightly bound disclina- 
tion dipole, leading to a defect density 
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As discussed elsewhere [l|, 0, 9 the number of disclina- 
tions, Ni, is determined by the Eulcr characteristic x of 
the background Ni = 6\, thus giving 12 for the sphere. 
We also note that the Burgers vector b is perpendicular to 
the dipole direction defined by the vector connecting the 
two disclinations forming the dislocation. The elastic en- 
ergy, Eq. |4]), includes now a disclination-dislocation en- 
ergy Eud and a dislocation-dislocation energy Edd given 
by 



E 



Dd 



r r Nl i 



i=i 



JV 2 



3 = 1 

2 Ni N 2 



YR 2 



i=l ] = 1 



and 

Edd 



N 2 

= Y f I ^x^yjS^^x-x,) J- 




FIG. 4: (Color online) Schematic diagram of a single scar 
aligned along a geodesic meridian on the two-sphere. The 
red arrows indicate the associated Burgers vector for each 
dislocation. 
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In spherical coordinates, we have b l a e a pdp = bgd^ — 
We assume that both components of the angu- 
lar Burgers vector, b\ and b^, are approximately \b t \/R, 

where \b l \ is taken to be the lattice spacing a. An explicit 
expression for the energy of an arbitrary dislocation dis- 
tribution interacting with N± disclinations is provided by 
combining Eqs. {SJ and (JOj) . 

Let us consider geodesically straight scars, symmet- 
ric about their midpoint and aligned along the fixed- cf> 
meridian, as shown in Fig. 2] With this choice Eq. ^ 
gives fa — <f>j and cos (3 = cos(9i — 9j). The Burgers vec- 
tors ¥, are orthogonal to the disclination dipole and 

symmetry implies that J^- b l = 0. 

Since we shall only consider glide motion for which 
dislocations move in the <f> direction, we may set b l e = 0. 
The elastic D-d interaction, Eq. ([8]), then reduces to 
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with 
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The dislocation-dislocation interaction, Eq.©, becomes 

N 2 N 2 

E dd = Y,Y,£dd{Pa) , (12) 
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FIG. 5: (Color online) The Disclination-dislocation (D-d) and 
dislocation-dislocation (d-d) interaction energies as a function 
of defect separation ft The inset is a blowup of the short- 
distance region of the plot. 
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FIG. 6: (Color online) The disclination-dislocation (D-d) in- 
teraction energy versus angular separation (3 for a variety of 
system sizes. 
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From now on we confine ourselves to the interaction be- 
tween defects in a single scar and ignore the effects of the 
neighboring scars. In Fig. [5] we plot Zr> d (solid line) and 
£dd (dotted line) as a function of the angular separation 
ft in units otYa 2 . Note that the disclination-dislocation 
interaction is attractive (for sufficiently short angular dis- 
tance) while the dislocation-dislocation interaction is re- 
pulsive. The formation of grain boundary scars may now 
be understood as arising from the competition between 
the attractive binding of a dislocation to an excess discli- 
nation (D-d interaction) and the mutual repulsion be- 
tween dislocations (d-d interaction) . Figure [5] shows the 
D-d interaction energy as a function of angular distance 
for a variety of system sizes. The functional dependence 
of £Dd/(Ya 2 ) on R/a given by Eq. (fTTj) implies that the 
strength of the short-distance attraction increases with 
system size. As a result the strong D-d attraction for 
large systems leads to more excess dislocations within 
a scar to stabilize geometric frustration. Note that the 
crossover from an attractive to a repulsive interaction oc- 
curs at a universal value of the order of 1 rad, consistent 
with the predictions of Refs. 0, Q . 

The potential energy E s of a single scar can be now 
obtained directly from the finite sum of pair interactions 
between defects: 



(14) 
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where ft is the angular distance of the ith dislocation 
from the center of the scar. 



We now compute the elastic stiffness [the spring con- 
stant in Eq.{T])] of dislocations within a scar. For that 
purpose, we consider small fluctuations of dislocations. 
Let us consider now small fluctuations of dislocations 
around their equilibrium positions in a scar with a fixed 
central disclination. The geodesic displacement of the ith 
dislocations will be denoted by Sj. We assume that the 
ith dislocation glides along the direction defined by the 
geodesic that starts at a point 9 di and forms an angle of 
ir/2 with respect to the arc connecting 9 di to the other 
defect location under consideration (the disclination 9o 
or the jth dislocations 9 dj ). The deformed geodesic arc 
distances [3 will then be 
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where 9jj, 9 di , and 9 d are the initial locations of the 
disclination and the ith and jth dislocations, respectively, 
and Si,Sj <C 0. For an arbitrary scar along a meridian, 
the relation between (3 and [3 is given by 

cos ft = cos ft cos(si/R) , 
cosftj = cos ftj cos[(si — Sj)/R] . (16) 

Expanding to second order gives 
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The first derivatives are easily seen to vanish, confirming 
that the initial configuration (linear and central symmet- 
ric) is a local minima: 
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The second derivatives are given by 
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Equation (|17p can then be written in terms of effective 
spring as 

AE S = \ £ K lS \ + \ Kit ( Si - S] f , (20) 
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with /3j and /Jy determined by the initial configuration 
to be 



We note that the deformation energy in Eq. ([20]) 
contains nondiagonal terms induced by the dislocation- 
dislocation interactions. The two stiffness coefficients Ki 
and Kij result from D-d attractions (£dd) and d-d re- 
pulsions (£dd), respectively, which implies that Ki > 
and Kij < for sufficiently short angular distance — i.e., 
(a/R) < P < 1 rad. Summing up we may write the 
energy shift as a general quadratic polynomial 



AE S = - ^2 KijSiSj , 



(24) 



with ICij given by 




(25) 



For a pinned, small-angle grain boundary in flat space, 
the restoring force to shear stress has been obtained in 
Ref. [22| , where it results from dislocation-dislocation in- 
teractions alone. The presence of disclination-dislocation 
interactions is a special feature of the two-dimensional 
curvature of the crystal. The eigenvalues ki of the ma- 
trix /C, 



KV = hV 



(26) 



give the effective stiffness coefficients with negative val- 
ues, indicating that the associated dislocation will not 
bind to an equilibrium scar. 



III. COMPARISON TO EXPERIMENT 

We now use the formulas developed in the previous sec- 
tion to compute the elastic stiffness. The elastic stiffness 
depends on the particular configuration of dislocations. 
We compare our results with the experimental data given 
in Ref. [12]. The stiffness coefficients will depend on the 
detailed defect spacings in the ground-state configura- 
tions as well as the total number of excess dislocations in 
a scar. We assume that dislocations are equally spaced, 
as actually observed experimentally (see Fig. [5]), although 
theoretical calculations predict that the spacing should 
grow with distance from the center of the scar [2|, |j] . We 
take the first dislocation be a distance 2a from the cen- 
tral disclination with the remaining dislocations spaced 
a distance 3a apart, as shown in Fig. [71 
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We note that the expressions for Ki and Kij show sin- 
gularities at f3 = 0, 7T. Those singularities are not real, 
as the validity of the above expressions is limited to 
(3 > (a/R). Although expressions that correctly capture 
the P — v limit may be derived, they are not necessary 
for the subsequent analysis. 
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FIG. 7: (Color online) Dislocation spacings for a linear sym- 
metric scar as highlighted within the blue box in Fig. [T] 
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1.28 
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15.20 


2.52 


1.00 


0.44 






50 


18.19 


3.57 


0.98 


0.83 


0.07 



TABLE I: Comparison of the numerical spring constants (in 
units of KbT/o 2 ), for individual dislocations, with the exper- 
imental values [13 ] . 



For the numerical _spring constants, we use the experi- 
mentally measured two-dimensional Young's modu- 
lus Y = 167 ksT/ a 2 . Note that the units of ki and Y are 
all k B T/a 2 . 

Spring constants for typical experimental sizes are 
plotted in Fig. [8] We see clearly that the elastic stiffness 
falls quickly as a function of distance from the central 
disclination, in agreement with general experimental re- 
sults. At a quantitative level, the predicted values for 
the elastic stiffness, in unit of ksT/a 2 , are compared to 
the experimental quoted results of Ref. [HJ in Table [J 
The continuum model predicts a very large value for the 
stiffness of the dislocation closest to the central disclina- 
tion, and indeed, this dislocation appeared immobile in 
the experiments in Ref. , and no elastic stiffness could 
be measured. The results for the next dislocations are in 
very good agreement, more so given that the experiment 
only contains a single scar realization, and the theoret- 
ical calculations ignore interactions among scars or the 
coupling of dislocations to the underlying lattice. 



IV. DYNAMICS OF VACANCIES AND 
INTERSTITIALS IN SPHERICAL CRYSTALS 
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We now analyze the dynamics of interstitials and va- 
cancies. In Ref. 16], three types of vacancies Vi a (crushed 
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FIG. 8: (Color online) Spring constants ki for each dislocation 
di, calculated for the configuration shown in Fig. [7] 
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FIG. 9: (Color online) We first consider an initial lattice 
with icosahedron symmetry, an (8, 3) icosadeltahedral lattice 
(stage 0). The I 2 interstitial is generated by adding one extra 
particle (stage 1), which evolves into a I3 interstitial (stage 
2). The curvature-driven unbinding of dislocations starts — 
the decay of an interstitial (stages 3 and 4). Individual dislo- 
cation glides towards the nearest isolated disclination (stage 
5 and 6). Each dislocation binds to a disclination to form 
three miniscars (stage 7). The results are obtained from the 
java applet developed in Ref. [23[, according to the procedure 
described in the Appendix. 



vacancy), SV (split vacancy), and V3 (threefold symmet- 
ric vacancy) were identified together with three inter- 
stitials ^3 (threefold symmetric interstitial), I2 (twofold 
symmetric interstitial), and l^a (disjoint twofold symmet- 
ric interstital). The I3 interstitial and V3 were found to 
be the most stable. Subsequent experiments [ItJ showed 
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FIG. 10: (Color online) We first consider an initial lattice 
with icosahedron symmetry, an (8, 3) icosadeltahedral lattice 
(stage 0). The SV vacancy is generated by subtracting one 
particle (stage 1), which evolves into a S3 vacancy (stage 2). 
The curvature-driven unbinding of dislocations starts — the 
decay of an vacancy (stages 3 and 4). Individual disloca- 
tion glides towards the nearest isolated disclination (stages 
5 and 6). Each dislocation binds to a disclination to form 
three miniscars (stage 7). The results are obtained from the 
java applet developed in Ref. [23I ] . according to the procedure 
described in the Appendix. 



that the different interstitials and vacancies exist as sta- 
ble defects, and their dynamics has recently been stud- 
ied [18j | . This situation is in contrast with spherical crys- 
tals where, to our knowledge, no interstitials or vacancies 
have been observed. 

In order to investigate vacancies and interstitials we 
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0.3231 
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0.4224 
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0.3231 


0.9420 


0.8003 


0.5594 


0.5594 
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0.2116 


1.0488 


0.8767 


0.7379 


0.7379 
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0.1222 


1.0739 


0.9601 


0.9420 


0.9420 



TABLE II: Angular distance between defects during the re- 
laxation of the interstitial shown in Fig. [9] 




Stages 

FIG. 11: (Color online) Continuum elastic energy from 
Eq. (fH)) for the configurations in Figs. and [TO] The ini- 
tial configuration is an I3 interstitial or a Vz vacancy, both of 
which are unstable to scar formation. The small difference in 
elastic energy for both I3 and V3 results from different final 
configurations (see Figs. [9] and I10[) . 

consider a system of 972 particles interacting with a 
Coulomb potential. In the initial configuration, there 
are only 12 disclination defects with the symmetry of the 
icosahedron, as shown in Fig. [9] (stage 0); this is an (8,3) 
icosadeltahedral configuration. We now force an I2 inter- 
stitial by adding a particle, as shown in Fig. [5] (stage 1). 
I2 evolves into I3, the bound complex of dislocations with 
zero net Burgers vector (stage 2). It is found that the 
^3 interstitial is unstable and starts to be ripped apart 
into three dislocations (stages 3 and 4) and eventually 
becomes three separate dislocations which each glide to- 
ward a fivefold disclination (stage 5 and 6). They quickly 
form a miniscar (a 5-7-5 grain boundary) at each of the 
vertices 5s by joining the nearest disclinations (stage 7). 
Snapshots of the dynamical sequence discussed above are 
shown in Fig. [9] 

A similar analysis may be done for vacancies. By 
subtracting a particle from the initial icosadeltahedral 
configuration (8,3) in Fig. \TU[ the lattice develops the 
structurally unstable SV vacancy (stage 1), which sub- 
sequently evolves into V3 (stage 2). Due to the energetic 
instability, V3 eventually forms three scars via curvature 
driven unbinding (stages 3-7), similar to the interstitial. 

Similar results also follow by considering 972 particles 
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interacting with a general potential 1 /r 7 (the generalized 
Thomson problem) , thus showing that the instabilities of 
vacancies and interstitials are a universal feature of two- 
dimensional spherical crystals. 

The instabilities of vacancies and interstitials are pre- 
dicted from the continuum elastic model described in this 
paper. To apply Eq. (|14p we first need to estimate the an- 
gular distances between defects at each stage. An (8, 3) 
icosadeltahedral lattice, with M = 972 particles, corre- 
sponds to system size R/a » 8.2 [using M « ^j=(-^) 2 ]. 

The relevant angular distances between defects can then 
all be calculated by simple counting together with spher- 
ical trigonometry The results so obtained are shown in 
Table [III Taking the orientation of the Burgers vectors 
appropriately into account we may then compute the to- 
tal interaction energy at each stage using Eq. (fT4f . The 
evolution of the total energy is shown in Fig.[TT] The en- 
ergy monotonically decreases until the final scarred state 
is reached. 



V. CONCLUSIONS 

In this paper we have studied the dynamics of point 
defects in two-dimensional spherical crystals. Our re- 
sults provide explicit predictions for the dislocation elas- 
tic stiffness that compare quite favorably with experi- 
mental data. We also analyzed the dynamics of intersti- 
tials and vacancies and found that the effects of curvature 
are quite dramatic, as defects that are stable in flat space 
become unstable in curved space. 

A number of issues raised in this paper will require fur- 
ther investigations. The process of generating vacancies 
and interstitials can be repeated indefinitely It should 
be expected that in this way, we can grow longer scars, 
whose length should saturate at some point. What struc- 
tures follow will be the subject of further investigations. 
Although our presentation has focused on the sphere, the 
results should apply equally to other geometries. It is 
expected that in arbitrary geometries, vacancies or inter- 
stitials should become unstable to the formation of scars 
nucleated by existing disclinations. 

Detailed experimental verification of our results could 
be achieved from experiments of colloids absorbed on oil- 
water interface as in Refs. [ToL H^ |. Using holographi- 
cal optical tweezers [ijj applied to a spherical crystal, 
it should be possible to remove one colloid, thus creat- 



ing a vacancy, which according to the results in our pa- 
per would become unstable and join existing scars, which 
could be visualized as described in Ref. [12j| . More rig- 
orous validations for the predictions in this paper can be 
accomplished by a more comprehensive analysis of exper- 
imental data such as the one presented in Ref. [12j |. 

In summary, the results presented in this paper show 
the dramatic effects of curvature in two-dimensional crys- 
tals. It is our expectation that this paper will motivate 
further experimental and computational work. 
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APPENDIX: ANALYSIS OF VACANCIES AND 
INTERSTITIALS USING THE JAVA APPLET 

The analysis of the stability of vacancies and intersti- 
tials (see Figs.|9]and HOl) has been obtained using the java 
applet available at Ref. 23]. In order to reproduce the 
results, we generate a (8,3) icosadeltahedral tessellation 
using the Construct(m, n) algorithm with m — 8 and n. 
Now add or remove a single particle to the lattice at the 
barycenter of a spherical triangle whose vertices are three 
nearest-neighbor five-fold disclinations by (shift + click) 
or (Ctrl + click). The self-interstitial (or self-vacancy) 
so formed is then relaxed by a standard relaxation al- 
gorithm. One immediately finds that a V% interstitial 
(or a SV vacancy) is structurally unstable. In a few 
time steps it morphs into a complex of dislocations with 
zero net Burgers vector — the most common structure ob- 
served is a set of three dislocations (Is or V 3 ) arranged 
in a hexagon. Removing a particle (or adding a particle) 
back restores the particle number to the original 972 and 
relaxing still leaves scars with total energy lower than 
the starting configuration with 12 isolated 5's. This es- 
tablishes that scars are definitely low-energy equilibrium 
states rather than artifacts of the relaxation algorithm. 
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